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Abstract 

We review results on the spherically symmetric, asymptotically flat 
Einstein- Vlasov system. We focus on a recent result where we found 
explicit conditions on the initial data which guarantee the formation 
of a black hole in the evolution. Among these data there are data 
such that the corresponding solutions exist globally in Schwarzschild 
coordinates. We put these results into a more general context, and we 
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include arguments which show that the spacetimes we obtain satisfy 
the weak cosmic censorship conjecture and contain a black hole in the 
sense of suitable mathematical definitions of these concepts which are 
available in the literature. 
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1 Some general remarks on gravitational collapse 

We start these notes with some general, informal, and in part historical 
remarks on concepts related to the phenomenon of relativistic gravitational 
collapse. 

Shortly after A. Einstein published his theory of general relativity 
|24l I25j . K. Schwarzschild showed that the following metric solves the cor- 
responding field equations in vacuum [39] : 

ds^ = - (^1 - df + (^1 - ~' dr^ + r\de^ + sin2 Od^^). (1.1) 

Here M > is a parameter, t € M is a time coordinate, the spacetime 
is spherically symmetric, and the polar angles and ip coordinatize the 
surfaces of constant t and r > 0. The latter are the orbits of SO (3) which 
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acts isometrically on this spacetime, and Airr'^ is the area of these surfaces. 
The part of this metric with r > 2M can be thought of as representing 
the gravitational field outside a static, spherically symmetric body of mass 
M and radius larger than 2M. With this interpretation in mind the radii 
r = 2M and r = 0, where the metric looks singular, lie within the matter 
where the metric does not apply anyway so that one need not worry about 
this singular behavior. However, using the new time coordinate 

i:=t + r* where r* := r + 2M ln(r - 2M), (1.2) 

the metric takes the form 

which extends smoothly through r = 2M, i.e., r = 2M is only a coordinate 
singularity in On the other hand, the so-called Kretschmann scalar 

K := R'^^^'Rap^s (1.4) 

derived from the Riemann curvature tensor R^^^^ blows up like M'^/r^ as 
r ^ 0; Greek indices are running from to 3 and are summed over if they 
appear as both lower and upper indices in the same expression. Since K is 
a scalar quantity, its values do not change under a change of coordinates. 
The singularity at r = is a genuine feature of the Schwarzschild spacetime 
which cannot be cured away by a more judicious choice of coordinates, it is 
a spacetime singularity where the structure of spacetime itself breaks down. 
Although the surface r = 2M is not singular there is something special about 
it, as can be seen from ()1.3p . The line r = 2M, 6 = Tr/2, = represents 
the world line of a massless particle (a photon) moving radially outward, 
away from the origin. Since no material particle can move faster than light, 
particles and photons can only pass through the surface r = 2M inward, but 
can never leave the region r < 2M once they are inside. Such a surface was 
later termed an event horizon. One should keep in mind from the above that 
by switching to different coordinates one was able to extend the spacetime 
beyond the region covered by Schwarzschild coordinates, even though in the 
latter coordinates the metric blows up at the boundary r = 2M. 

While the part with r > 2M of the Schwarzschild metric was successfully 
used to explain for example the perihelion precession of Mercury in the solar 
system, it was argued that no conceivable physical process could compress 
an amount of matter so much that its mass M would all be contained inside 
the region r < 2M, and so the irritating behavior at the surface r = 2M and 
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also the break down of the geometry of spacetime at r = was discarded as 
unphysical. 

But in 1939 J. R. Oppenheimer and H. Snyder [32] constructed a semi- 
exphcit, time dependent solution of the Einstein field equations where a 
homogeneous spherically symmetric ball of dust, i.e., of a fluid with pres- 
sure zero, collapses until all the mass is within the region r < 2M, and it 
continues to collapse until the scalar curvature of spacetime blows up at the 
center of symmetry. Although this proved that an event horizon can evolve 
out of completely regular initial data, several decades passed before such 
structures became accepted as potentially relevant from the physics point 
of view and J. A. Wheeler coined the name black hole for them. 

Today, there are many astronomical observations for which the currently 
best explanation involves a black hole. For example, black holes of the 
order of 10^-10^ solar masses are believed to reside in the centers of many 
galaxies, including the Milky Way. In spite of the increasing relevance of 
black holes as real astrophysical objects many important basic questions 
about gravitational collapse are still open. The most prominent of these is 
the cosmic censorship conjecture. 

To see the issue here, we notice that in the Oppenheimer-Snyder example 
the spacetime singularity which forms at r = is hidden behind the event 
horizon so that it cannot be seen or in any other way be experienced by 
observers outside the event horizon. The same is true for the spacetime sin- 
gularity in the Schwarzschild spacetime. In the 1960s R. Penrose [33] proved 
that a spacetime singularity forms in the gravitational collapse of a not nec- 
essarily symmetric star made up of "reasonable" matter, i.e., matter which 
satisfies the strong energy condition, provided a closed trapped surface forms. 
This is a closed, two-dimensional, spacelike surface with the property that 
the null geodesies, i.e., the light rays, which start perpendicular to the sur- 
face decrease its surface area both when followed inward and outward from 
the surface; for a precise definition cf. |51[ 9.5]. The surfaces of constant t 
and r with r < 2M in (II. 3p are trapped, and we note that Schwarzschild co- 
ordinates cannot cover regions of spacetime which contain trapped surfaces. 
Since trapped surfaces are stable under small perturbations of the spacetime, 
Penrose's result showed that the formation of spacetime singularities is not 
restricted to spherically symmetric, especially constructed or isolated exam- 
ples but is a genuine, stable feature of spacetimes. However, his result gives 
little information about the geometric structure of a spacetime with such a 
singularity, in particular it does not provide an event horizon surrounding 
the singularity. The existence of a naked singularity which by definition is 
not hidden behind an event horizon would violate predictability as it would 
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not be possible to predict from initial data what an observer would see if he 
could observe a singularity. Hence Penrose formulated the cosmic censorship 
conjecture which demands that any singularity arising in the gravitational 
collapse of generic regular initial data is hidden behind an event horizon; 
exceptional data leading to naked singularities are required to form a "null 
set" in some suitable sense. The above is an informal statement of the so- 
called weak cosmic censorship conjecture [51j 12.1]. It would in particular 
guarantee that predictability holds at least in the region outside the event 
horizon. In the strong version no observer is allowed to observe a singularity. 
For a mathematical definition and discussion of the weak and strong cosmic 
censorship conjectures we refer to |20j . 

An important example where naked singularities do form for a null 
set of data while cosmic censorship holds for generic data, is the spher- 
ically symmetric Einstein-scalar field system which was investigated by 
D. Christodoulou, cf. [Ml [19]. A massless scalar field or dust as employed by 
Oppenheimer and Snyder and later also by Christodoulou are but two 
possibilities for modeling matter in gravitational collapse. In the present 
notes we discuss results where a collisionless gas as described by the Vlasov 
equation is used as matter model, a model which we consider particularly 
suitable for this purpose from a mathematics point of view and which is well 
motivated from an astrophysics point of view, cf. [llj. 

2 The Einstein- Vlasov system 

Consider a smooth spacetime manifold M equipped with a Lorentzian metric 
gaf3 with signature ( — h -|- -|-). The Einstein equations read 

= SnTaiB, (2.1) 

where G^/j is the Einstein tensor, a non-linear second order differential ex- 
pression in the metric Qa/s, and T^^ is the energy-momentum tensor given 
by the matter content (or other fields) of the spacetime. To obtain a closed 
system, the field equations (|2.ip have to be supplemented by evolution equa- 
tion(s) for the matter and the definition of T^/s in terms of the matter and 
the metric. 

We choose as matter model a collisionless gas. In order to write down an 
evolution equation for the number density of the particles on phase space we 
recall that the world line of a single test particle on M obeys the geodesic 
equation 

i" = p", p" = -nyp^ (2.2) 
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where denotes general coordinates on M, are the corresponding canon- 
ical momenta, F^^ are the Christoffel symbols induced by the metric 5^^, 
and the dot indicates differentiation with respect to an affine parameter, 
i.e., with respect to proper time along the world line of the particle. We 
assume that all the particles have the same rest mass, normalized to unity, 
and move forward in time. Hence, their number density / is a non-negative 
function supported on the mass shell 



PM := [g^ppy^ = -1, > 0} 



a submanifold of the tangent bundle TM of the spacetime manifold M which 
is invariant under the geodesic flow. Letting Latin indices range from 1 to 
3 we use coordinates {t^x"") with zero shift which implies that g^a = 0. On 
the mass shell PM the variable then becomes a function of the remaining 
variables (t , , p^) : 



p 



P ^ y^Vl+5a6pV- (2.3) 



Since the particles move like test particles in the given metric, their number 
density / = f{t,x°',p^) is constant along the geodesies, and it satisfies the 
Vlasov equation 

dtf + ^ d,af - 1 r^y^^ dpaf = 0. (2.4) 

The energy-momentum tensor is given by 

rr f f\ 11/2 dp^dp^dp^ 
J -po 

where \g\ denotes the modulus of the determinant of the metric, and indices 
are raised and lowered using the metric, i.e., pa = gapP^ ■ The system 
(j2.ip . (j2.4p . (j2.5p is the Einstein- Vlasov system in general coordinates. An 
introduction to relativistic kinetic theory and the Einstein- Vlasov system 
can be found in [H [37]. The Vlasov equation is widely used as a matter 
model in astrophysics, to describe galaxies or globular clusters [H]. Such 
systems are usually dealt with as isolated systems in an otherwise empty 
universe which in our context means that the spacetime is asymptotically 
flat. 

Let us for a moment consider a distribution function of the form 

f{t, x'^,p') = -uo \g\-'/^p{t, X-) 5{p' - u\t, X-)), (2.6) 
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where p = p{t, x"") is a scalar function on spacetime, 5 is the Dirac 5- 
distribution, and = takes values in the mass shell so that uq 

is determined by v!^ , cf. (j2.3p . Then the macroscopic quantities p and 
together with the metric satisfy the Einstein-Euler system for a perfect fluid 
with pressure zero, a matter model referred to as dust in the first section and 
used by Oppenheimer and Snyder; p is the mass-energy density of the dust- 
fluid and its four-velocity. It should be stressed that although formally 
the Einstein-dust system can be viewed as a special case of the Einstein- 
Vlasov system we use the term Vlasov matter exclusively for genuine (and 
usually smooth) distribution functions on the mass shell PM. 

Before we proceed a few advantageous features of the Einstein- Vlasov 
system are worth to be pointed out. Firstly, if the metric and therefore 
the Christoffel symbols are given, the evolution equation for the matter, 
i.e., the Vlasov equation, does not produce any singularities by itself; it 
is — for a given metric — indeed a linear first order conservation law which 
can be solved by the method of characteristics. This situation is different if 
the matter is described as a fluid, and singularities induced by the matter 
model can prevent one from analyzing the formation of event horizons and 
true spacetime singularities. Secondly, in the Newtonian limit the Einstein- 
Vlasov system turns into the Vlasov-Poisson system \AQ\ H6] for which there 
is a global existence result for general, regular data, cf. [301 EH EH] • Hence 
any singularity in the solutions of the Einstein- Vlasov system should have 
its origin in some relativistic effect; for a fluid this is again different. Beside 
these mathematical properties the Vlasov equation also has a clear physical 
interpretation and motivation, as pointed out above. 

The questions raised in the previous section are at present out of reach 
of a rigorous mathematical treatment, unless simplifying symmetry assump- 
tions are made. Hence we will consider the Einstein- Vlasov system under 
the assumption of spherical symmetry. Notice that the investigations which 
we cited above and where dust or a massless scalar field are used as mat- 
ter models employ the same symmetry assumption. We use Schwarzschild 
coordinates (t, r, 6, ip) and write the metric in the form 

as to the range and meaning of these coordinates we refer to the previous 
section, cf. (jl.ip . Asymptotic flatness means that the metric quantities A 
and /i have to satisfy the boundary conditions 

lim X{t,r) = lim fi{t,r) = (2.7) 

r— >oo r— >oo 
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so that as r — >■ oo the metric approaches the fiat Minicowsici metric, written 
in polar coordinates. In addition we impose the boundary condition 

X{t,0) = (2.8) 

in order to guarantee a regular center. Since polar coordinates sometimes in- 
duce artificial coordinate singularities at r = 0, it is convenient to introduce 
the corresponding Cartesian coordinates 

X = {x^,x'^,x^) = r {sin COS if, sin sin if, COS 9). 

li p = ip^,p^,p^) denotes the corresponding canonical momenta, then 

where |pp = (p^)"^ + (p^)^ + (p^)^ and x ■ p = x^p^ + x^p^ + x^p^. Since 
this quantity appears in the formula for the components of the energy- 
momentum tensor which in turn appear as source terms in the field equa- 
tions, it is preferable to use the non-canonical momentum variables 

^;«=p« + (e^-l)^^, a = 1,2,3. 
r r 

In these variables, 

Po = -e^Vl + bl^ 
and / is spherically symmetric iff 

f{t,x,v) = f{t,Ax,Av), x,^;GM^ AeS0(3). 

The spherically symmetric, asymptotically flat Einstein- Vlasov system takes 
the following form: 

dtf + e^-^^^== ■ d^f - (dtX^ + e^^-^fi^/TTW) --9^ = 0, 



1 -f- \vY' \ r J r 

(2.9) 

e-'^^(2rdr\ -!) + ! = 87rrV, (2.10) 
e-^^{2rdrfi + 1) - 1 = STrr^p, (2.11) 
dtX = -Airre^+f^j, (2.12) 



(^d^H + {drii - drX){drfi + ^)^ - e'^^^ {dfx + dtX {dtX - dtp.)) = Airq, 

(2.13) 
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where 



p{t,r) 



p{t,x) = / -v/l + \v\'^f{t,x,v) dv 



(2.14) 



p{t,r) 



f / X ' V \ ^ 

p{t,x) = / [—^) f(t,x,v) 



dv 



(2.15) 



j{t,r) 



q{t,r) 




dv 



(2.16) 



(2.17) 



For a detailed derivation of these equations we refer to [37]. It should be 
noted that in this formulation no raising and lowering of indices using the 
metric appears anywhere. It is a completely explicit system of PDEs where 
x,v G M^, X ■ V denotes the Euclidean scalar product, and = v ■ v. 

For the spherically symmetric Einstein-dust equations Christodoulou |13j 
showed that cosmic censorship is violated. Indeed, not even a suitable small- 
ness condition on the initial data prevents the formation of naked singular- 
ities for dust. As we shall see in the next section this is different for the 
Vlasov matter model, provided we have a genuine, smooth distribution func- 
tion with respect to x and v. 

The goal for the spherically symmetric, asymptotically flat Einstein- 
Vlasov system is to show that for all (or at least for all generic) regular 
initial data the corresponding solution is either global in the sense that 
the spacetime is singularity-free or the solution undergoes a gravitational 
collapse in which a spacetime singularity forms which is hidden behind an 
event horizon. Of course the ultimate goal would be to prove this for general, 
not necessarily symmetric data, but let's be modest for now. So far all 
analytical and numerical results support the conjecture that the above is 
indeed true and that in particular the spherically symmetric, asymptotically 
flat Einstein- Vlasov system satisfies the weak cosmic censorship conjecture. 
An existence result for singularity-free solutions for restricted, small data 
has been known for some time and is reviewed in the next section. A class 
of data which lead to gravitational collapse as described above has been 
established more recently, and this result is discussed in Section [H 

3 Local and global existence results 

In this section we review a number of results from the literature for the 
spherically symmetric, asymptotically flat Einstein- Vlasov system. They 
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serve as either background or counterparts to the results on gravitational 
collapse which we state and discuss in the next section. 



3.1 Local existence and continuation 

Due to our choice of w as a non-canonical momentum variable the system 
(j2.9P " (j2.17p has the following nice feature. If a distribution function / is 
given, then the source terms /o, p, j, and q can be computed from it without 
reference to the metric. Given p and observing the boundary condition (j2.8p 
the field equation (j2.10p can be integrated to yield 

o\ 2m 

= 1 , (3.1 

r 

where the quasi-local mass m is given by 

IT 

m{t,r) := 4-7T / p{t,ri) ■q'^drj. (3-2) 



Given p and A the field equation (12. lip together with the boundary condition 
()2.7p determines p: 



oo 



p{t, r) = exp ( - / e2^(*''') ( + 47Tr,p{t, t^)) drj) . (3.3) 



Finally, if both A and ^ are given and sufficiently regular, then / is deter- 
mined from its initial data by the method of characteristics: 

f{t,x,v) = f{{X,V){0,t,x,v)) (3.4) 

where {X,V){s,t, x,v) is the solution of the characteristic system 

X = e^-^ -^^== ,v = -(dtX^ + e'^-'^drfi^l + W) - (3-5) 
■s/l + V r J r 

o 

of the Vlasov equation (j2.9p satisfying {X,V){t,t,x,v) = {x,v), and / is 
the prescribed data at time t = 0. Notice that the characteristics are now 
parameterized by coordinate time instead of proper time as in (j2.2p . 

The iterative scheme indicated above can be used to prove a local exis- 
tence and uniqueness theorem together with an extension criterion. 

Theorem 3.1 Let f£ C^(M^) be non-negative, spherically symmetric, and 
such that for the induced quasi-local mass, 

2rh{r)/r < 1, r > 0. (3.6) 
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Then there exists a unique regular solution f of the asymptotically flat, 
spherically symmetric Einstein- Vlasov system with /(O) = f on a maximal 
interval of existence [0,T[ with T > 0. // 



then T = oo. 

This result was first established in [39], see also [37]. Some comments are in 
order. 

Functions in C^{M.^) are by definition continuously differentiable and 
compactly supported. A solution is called regular if the derivatives which 
appear in the system exist in the classical sense and are continuous; for the 
precise definition we refer to [37J. 

The restriction (|3.6p on the initial data is necessary in view of (|3.ip 
because that equation defines A only as long as 2m{t, r)/r < 1 so we have to 
require this for the initial data. This is related to the fact that, as noted in 
the first section, Schwarzschild coordinates cannot cover regions of spacetime 
which contain trapped surfaces. 

A local existence and uniqueness result for the Einstein- Vlasov system 
without a symmetry assumption was established by Y. Choquet-Bruhat }12j . 
However, in order to extend a local solution to a global one based on the 
latter result one would have to control high order Sobolev norms of the 
solution. The extension criterion provided in Theorem 13.11 is much less 
demanding and forms a more convenient starting point for investigating 
global properties of solutions. 

The major simplification of the system due to the symmetry assumption 
is the fact that for given source terms the metric is completely determined by 
the two constraint equations ()2.10p and (I2.11|) which are ordinary differential 
equations in r. The metric has no independent degrees of freedom, and this 
rules out gravitational radiation, much like the Vlasov-Maxwell system being 
reduced to the Vlasov- Poisson system if spherical symmetry is assumed. But 
if A is defined by (|3.ip . then it becomes rather unpleasant to control dtX 
which appears in the Vlasov equation and to make sure that the condition 
2m{t,r)/r < 1 is preserved. In [37] this was resolved by considering first a 
modified system where (|2.13p was left out and dtX was replaced by a quantity 
A defined by (j2.12p . Then one can show a posteriori that indeed dfX = A 
and that p.l3p holds as well. Eqn. (I2.12p can then be used to control dtX 
and therefore A in terms of 



sup< \v\ I {t,x,v) £ supp /, 0<t<r><oo 



P{t) := sup< 1^1 I {s,x,v) G supp/, < s <t 
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notice that by (j3.4p the function f{t) is bounded so that the source terms 
p,p,j, q are bounded by powers of P. 

Eqn. (j2.13p also has its role to play. To establish the convergence of the 
iterative scheme one needs uniform bounds on certain second order deriva- 
tives of the metric coefficients. But since in (j2.10p - (|2.12|) the first order 
derivatives depend in a pointwise way on the source terms, it seems that 
one needs to control derivatives of these source terms. In the corresponding 
results for the Vlasov-Poisson or Vlasov-Maxwell systems one exploits the 
fact that the field quantities depend on the source terms through spatial 
or spacetime integrals; the corresponding field equations are smoothing, cf. 
[lOl [26| [38] . Here the field equation (j2.13p provides this smoothing effect, 
because it turns out that in order to control the derivatives of the charac- 
teristic flow with respect to the initial data only a certain combination of 
second order derivatives of the metric coefficients is needed, and this combi- 
nation is precisely the one which appears in ()2.13p and is therefore controlled 
by q- 

In the context of these arguments and also in what follows some further 
information is useful. By (I2.10p and ()2.1ip . 

dr^i + drX > 0, 

and together with the boundary condition (|2.7p and (|3.ip this implies that 



Solutions of the Einstein- Vlasov system satisfy the following conservation 
laws. The Vlasov equation implies that the number of particles 



is conserved. More importantly, if we observe that j is for each fixed time t 
compactly supported so that by (j2.12p . dtX{t,r) = for r sufficiently large, 
we conclude from (13. 1[) that the ADM mass 



is conserved. 

3.2 Global existence for small data 

A natural question when investigating a nonlinear system of PDEs is whether 
small initial data lead to global solutions which disperse. As mentioned 



/U-A</x-hA<0. 



(3.7) 




(3.8) 




(3.9) 
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above, no such smallness condition can be formulated for the at least for- 
mally closely related Einstein-dust system, cf. [13] , but for true Vlasov mat- 
ter such a global result for small data does exist. 

Theorem 3.2 For all R > there exists e > such that if f is a maximal 
solution of the asymptotically flat, spherically symmetric Einstein- Vlasov 

o 

system with /(O) = / satisfying 



and 



f{x,v) = for \x\ + \v\ > R 



then the solution exists globally in t. Moreover, the solution disperses in the 
sense that 



\\pmoo<c{i + \t\) 

and the spacetime is geodesically complete. 



t G 



This result was first proven in [39], see also \37\ I47j . Similar results were 
known both for the Vlasov-Poisson and the Vlasov-Maxwell system, cf. f9l 
[27tl38j. The basic dispersive mechanism in the system can be seen as follows. 
By (|3.4p and the change of variables formula. 



p{t,x)= / y^l+WfiiX,V)iO,t,X,v))dv 



Vl + \v\^f{X,V{0,t,x,v)) 



det 



fdX{0,t,x,' 



dv 



dX. 



Now in the flat Minkowski case, A = /i = 0, X(0, t,x,v) = x — tv/ y^l + 
and the determinant above grows like t^. And if the field terms satisfy 
suitable decay conditions this turns out to remain correct. The decay of p 
and the other source terms in turn implies decay for the field terms, and a 
bootstrap argument gives the result. Here it again becomes important that 
a certain combination of second order derivatives of the metric coefficients 
can be controlled via the source terms by the field equation p.l3p . 

The above argument proves that the solution is global with respect to 
the chosen time coordinate. However, this does not automatically imply 
that the maximal Cauchy development of the corresponding initial data is 
singularity-free. By definition, a spacetime contains a singularity if there 
exists a timelike or null geodesic, i.e., a world line of a particle or a photon, 
whose maximally extended domain of affine parameter, i.e., proper time in 
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case of a particle, is not the whole real line. The geometrically invariant, 
coordinate-free characterization of a singularity-free or global spacetime is 
therefore that all timelike and null geodesies exist on the whole real line. In 
the course of the proof of the above bootstrap argument one obtains sufficient 
decay information on the Christoffel symbols to conclude that indeed all 
maximally extended geodesies in this spacetime are complete, i.e, exist on 
the whole real line. 

A mechanism which also leads to global existence but is different from 
dispersion induced by small data is that initially all the particles move out- 
ward sufficiently fast to prevent re-collapse, cf. [5j. For a collisionless gas of 
massless particles a result analogous to Theorem 13.21 has been shown in [22j . 

3.3 Possible break down must occur at the center first 

The next question is whether large data lead to a break down of the so- 
lution. It is important to realize that the fact that — as shown in the next 
section — certain data do lead to gravitational collapse and the formation of 
black holes does certainly not rule out the possibility that in Schwarzschild 
coordinates all solutions of the spherically symmetric, asymptotically flat 
Einstein- Vlasov system are global. If the weak cosmic censorship conjecture 
holds for this system, then, physically speaking, an observer who is very far 
away (ideally at spatial infinity) should not be able to observe the singu- 
larity so that for him the universe should look singularity-free, and in this 
context one should note that Schwarzschild time asymptotically coincides 
with the proper time measured by such an observer. So far, all numerical 
simulations and analytical results support the conjecture that all solutions 
are global in Schwarzschild time. In all cases where a gravitational collapse 
was simulated numerically the solution did not seem to break down in finite 
Schwarzschild time and an event horizon did form, cf. [71 [3T| . In the 
next section we discuss a recent analytical result which shows that at least 
for a certain class of data solutions do again exist globally in Schwarzschild 
time, but they do undergo a gravitational collapse and form a black hole. In 
the context of that result it is important to know that if a solution should 
develop a singularity in finite Schwarzschild time, then this must happen at 
the center first. 

Theorem 3.3 Let f be a regular solution of the spherically symmetric, 
asymptotically flat Einstein- Vlasov system on a time interval [0, r[. Sup- 
pose that there exists an open neighborhood U of the point (T, 0) such that 

sup{[t;[ I {t,x,v) G supp/ n ([/ X M^)} < oo. (3.10) 
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Then f extends to a regular solution on [0,T'[ for some T' > T. 



This theorem was estabhshed in |43j. In order to explain the mechanism 
behind this we introduce coordinates on the mass shell which are adapted 
to the symmetry: 



X ■ V 

\x\, w = , L = \x X V 

r 



2 



Spherical symmetry of / implies that by abuse of notation, 

f{t,x,v) = f{t,r,w,L). 

One should think of w as the non-canonical radial momentum and L as 
the modulus of angular momentum squared of a particle. Due to spherical 
symmetry L is conserved along characteristics, and the characteristic system 
written in (r, w, L) takes the form 

r = e'^"^|, (3.11) 

w = -(^dtXw + ef'~^drfiE-ef'-^-^^ , (3.12) 

L = 0, (3.13) 

where 



E = E{r, w, L) := y^l + w'^ + L/r^ = eV- 
The source terms take the form 

p[t,r) = — / Ef{t,r,w,L)dLdw, (3.14) 

J~oo Jo 

= / -^f{t,r,w,L)dLdw, (3.15) 

j{t^r) = — / w f{t,r,w,L)dLdw- (3.16) 
f J-oo Jo 

the quantity q is not needed in what follows. Since on the support of / 
the quantity L is bounded initially, it remains bounded for all time. In 
order to establish the theorem one needs to control the increase in |t;| along 
characteristics which stay away from the center. Because of the relation 
= nP' -\- L/r"^ it is sufficient to control w along such a characteristic. If 
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we substitute the expressions for dtX and drU and for the source terms into 
([3l^ . we find that 



S-g- 1 f{s,r,w,L)dLdw 

with obvious definition of E. Now L is bounded, r is bounded away from 
zero, m is bounded by M, and /i — A</i + A<Oby ()3.7p . Hence the term 
to focus on is the double integral. But 

WW — ii— = — = 

EE wE + wE 

_ w w'^il + L/r'^) - + L/r^) 

E wE + wE 

provided the denominator does not vanish. In the last step the term w'^nP' , 
which is the worst one concerning powers of w or w, canceled, and this is 
the crucial observation in order to establish Theorem 13.3^ since it leads to 
an estimate of the form 

P{t) < P(0) +C [ P(s) InP(s) ds, 
Jo 

cf. the continuation criterion in Theorem 13. 1[ 

The result above refers to the Einstein- Vlasov system written in 
Schwarzschild coordinates. The analogous result in so-called maximal- 
isotropic coordinates was established in [47J. In |23j the analogous result 
is established using coordinates which cover the whole spacetime, even in 
the presence of trapped surfaces. 

To appreciate this result it is instructive to think of a spherically sym- 
metric solution of the Einstein-dust system. In that system all particles 
move radially, and particles at the same radius r have the same momen- 
tum, i.e., they remain on the same sphere which can grow or shrink as a 
whole, cf. (j2.6p . If two such spheres cross, the system experiences a so-called 
shell-crossing singularity which occurs at some positive radius. If one wants 
to study the formation of event horizons and spacetime singularities in the 
Einstein-dust system, one has to handle the problem that the solution may 
break down due to a shell-crossing singularity before the objects of interest 
have actually evolved. 



w = e"-^^ - e^^^ "~'^:'' E 
r-^E 



m 








oo 


POO 


-oo 


Jo 



WW 
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4 Gravitational collapse and formation of black 
holes 



As pointed out repeatedly above, Schwarzschild coordinates do not cover 
regions of spacetime which contain trapped surfaces. Since in gravitational 
collapse the latter typically appear before a singularity forms, one may well 
argue that in order to investigate gravitational collapse and the formation 
of black holes, one should better not use these coordinates to begin with. 
On the other hand, they do have advantages from an analysis point of view, 
and one may still hope to derive the desired information on gravitational 
collapse from the asymptotic behavior of the solution for large Schwarzschild 
time. But for large data which possibly lead to gravitational collapse there 
is so far no global existence result in Schwarzschild time. In order to by-pass 
this difficulty we study the solutions in a coordinate region which avoids the 
center, on which we can consistently formulate and study the corresponding 
Cauchy problem, and which is large enough to allow us to conclude the for- 
mation of black holes. That the behavior of the solution on this coordinate 
region indeed implies the weak cosmic censorship conjecture and the forma- 
tion of a black hole is shown in Section [4.31 where we rely on a mathematical 
formulation of the corresponding concepts given in [20j . 



4.1 The system on an outer domain 

We study the solutions to the spherically symmetric, asymptotically flat 
Einstein- Vlasov system (l2.9p ~( f2?T7l) on the exterior region 

D:={(t,r) € [0,oo[2|r>7+(t)}, (4.1) 

where 7''' is an outgoing radial null geodesic originating from some r = rg > 
0, i.e., ^ 

^(s) = e{/^-A){.,7+(s))^ ^+(0) ^ (4_2) 
ds 

This equation is obtained from (j3.1ip by replacing the 1 in the definition of 
E by 0, i.e., the particle under consideration has rest mass like a photon 
should, and by putting L = 0, i.e., the photon moves radially. That the 
photon is outgoing, i.e., moving away from the center, means that w > 0, 
and hence w/E = 1. In order to restrict the analysis of the system to the 
region D we have to find a replacement for the boundary condition (j2.8p . 
We prescribe the total ADM mass M > and redefine the quasi-local mass 
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by 

m{t,r) = M - 4tt I p{t,r]) rj'^dr], (4.3) 



while retaining the definition (|3.1|) for A. Clearly, a solution of the system as 
considered in the previous section, when restricted to D, is a solution of this 
modified system. Moreover, characteristics of the Vlasov equation can pass 
from the region D into the region {r < 7^(i)} but not the other way around 

o 

so that initial data / posed for r > completely determine the solution on 
the outer domain D. Such data need to satisfy the restrictions (13. 6|) and 

/■oo 

Mont = 4^ / A??) r]^dr] < M, (4.4) 

o . ° 

where p is induced by /. Then lim^^oo Tn{t,r) = M and < m < M. 
The crucial question is whether one can specify data such that 7"^ has the 
property that 

lim 7^(s) < 00. (4-5) 

s— >oo 

While this is not sufficient to conclude the formation of a black hole, it turns 
out to be the main step towards that goal. 

It is important to note that the behavior or even the nature of the matter 
in the region {r < "y~^{t)} is not going to be relevant in what follows. For 
example one can equally well think of Vlasov data being posed on {r > 0}, 
but only the data on {r > tq} need to be properly restricted in order to 
obtain the desired behavior of the solution on the outer domain D. What 
is essential is that there initially is and hence remains some matter in the 
region {r < 7^(t)} as guaranteed by the condition Mout < M. 

4.2 The main result — analysis in Schwarzschild coordinates 

The initial data / £ C^{M.^) for the outer matter should satisfy the condition 

o 

that on supp/, 

Ro<r <Ri, w<W- 
where we use the variables (r, w, L) introduced above, and 

< ro < i?o < ^1, < 0. 

In particular, all particles move inward initially. These parameters can be 
specified in such a way that the data are close to violating the condition 
2rh(r)/r < 1, that the particles continue to move inward on the outer domain 
D, and that (14. 5p holds. The main result is the following. 
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Theorem 4.1 There exists a class of regular initial data for the spherically 
symmetric Einstein- Vlasov system such that for such data the corresponding 
solution exists on D, and 

lim 7^(s) < oo. 

In addition the following holds: 

(a) Fort ^ oo no matter remains in the region {r > 2M}, more precisely, 
the solution is vacuum and the metric equals the Schwarzschild metric 
with mass M in the region 

t>0 andr>2M + ae~^K 

Here a, /? > depend only on the initial data class. 

(h) In the outer region D and for r < 2M , 

lim ^{t, r) = — oo, 

more precisely, 

, , / ae-P' V'^ 
n(t, r) < In — 37 

for all t > and 7^(t) < r < 2M + ae~^^. This implies that for 
c < 2M the timelike lines r = c are incomplete, i.e., they have finite 
proper length, and this length is uniformly bounded. 

(c) The radially outgoing null geodesic which does not escape to infinity 
and is furthest to the right with this property gets trapped precisely at 
the Schwarzschild radius r = 2M . More precisely, we define 

r* := sup{r > tq | the radially outgoing null geodesic 7 with 
7(0) = r satisfies lim 7(5) < 00}, 

s— ►oo 

and let 7* be the radially outgoing null geodesic with 7*(0) = r* . Then 

lim 7*(s) = 2M, 

s— ►oo 

and every radially outgoing null geodesic 7 with 7(0) > r* is future 
complete, i.e., it exists on [0,oo[ in an affine parameterization, and 
lims^oo7('S) = oo- 
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If a solution with mass M of the sphericahy symmetric, asymptotically flat 
Einstein- Vlasov system collapses to a black hole of mass M and if we co- 
ordinatize this solution by Schwarzschild coordinates, then the asymptotic 
behavior as i ^ oo of this solution should be as given by this theorem. In 
the next subsection we show that this behavior actually implies that a black 
hole forms in the sense of a suitable, coordinate-independent formulation of 
this concept. 

For a complete proof of Theorem 14.11 we refer to [6] . Here we want to 
highlight some central arguments. The first of these makes sure that the 
particles in the outer domain D keep moving inward in a controlled way. 
Since initially w < W- < for all particles, this remains true on some time 
interval. On this time interval and along any characteristic in supp/, 

d , ^ , 4tt^ rco roc 
ds r 



oo Jo 

,,m f 1 + L/r"^ 

— 1- , , ^ 

r'^\ E r^E J r^E 

< -e'^^ ^ — + ^tt; + e^" 




2 1 E r^E r^E' 



Differentiating ()3.ip w.r.t. t and using (f2T2]l leads to dtm = -A^Tr'^e^'-^j > 
on the time interval we consider. It follows that m{s,r) > m(0,r) = m{r). 
Thus as long as the characteristic remains in 



ds 

We require that 



^(e "^tt;) < I L — — m(r) — r m(r) ) . 

is r-^E \ r J 



3L 

< L < — m(r) + r rn(r), r G [ro, (4-6) 
r 



on supp/. Then the above estimate together with a bootstrap argument 
and the fact that > 1 show that 

w<( min e"^^"''")^ W- (4.7) 

\r£[Ro,Ri] ) 

on supp f n D. 

Before going further some comments on the condition (j4.6p are in order. 
By our general set-up some mass must be inside {r < tq} initially, and this 
mass is a lower bound on m in ()4.6p . The role of this mass is not to pull the 
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particles inward, but to keep them focused towards the center. The smaUer 
their angular momentum is, i.e., the better they are aimed straight towards 
the center, the smaher can the central mass be chosen initially. Notice that 
for spherically symmetric dust which is used as matter model in [32] and 
|13| all particles have angular momentum equal to zero. 

The second issue in the proof of Theorem 14.11 we want to touch upon is 
the limiting behavior of 7'''. The basic idea is to consider a radially ingoing 
null geodesic 7" which starts to the left of the outer matter and to the right 
of 7+, i.e., 

^(s) = _e(/^-^)(«'7-W), ro<n= 7-(0) < Rq. 
as 

Then initially and therefore as long as 7^ and 7^ do not intersect there is 
no matter in the region between the outgoing and the ingoing null geodesic. 
This fact can be used to estimate their relative velocity in such a way that 
in Schwarzschild time they actually never intersect. This proves the limiting 
behavior of 7'*' and furthermore shows that the matter which is initially in D 
stays strictly to the right of 7^ and therefore in D for all future Schwarzschild 
time. In order to estimate how far the two null geodesies can move at most 
we observe that by (13. 3p . 

Kt,r) <- J\^^(''^^^^^ds =: ii{t,r), 

and hence 

|7±| < e^. 

We wish to see that the right hand side becomes small, and to this end we 
observe that the following chain of estimates yields a lower bound for dtfi- 

/oo 

/•oo 

= 47r / r]{p + p){t,ri)e^'^^e^^ dri 

Jr 

/oo 
3r]{-j{t,7]))e''+^e'^^ dr] 

<-3Ridtfi{t,r). (4.8) 

In the first estimate we exploit the fact that not only is j negative, but by 
choosing \W- \ large the source terms p and p can be estimated by a suitable 
multiple of —j. 
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If one puts all the details which are left out here together, one obtains 
a list of conditions on the initial data which make sure that these estimates 
hold true for all future Schwarzschild time. The resulting conditions are the 
following: 

o 

supp/C [Rq, Ri]x] — oo, Wo] X [0, oo[ 

with 

<ro <ri = 2M < Ro= ^ < Ri, W- < 0, 

0<niro)<M, (4.9) 



and either 



Ri-ri< 



or 



ri - rp 
6 ' 

^2 



ri (1 r-i - ro 

1 < mm < — 



Ri [6 36RiM 24i?i 

Since ri = 2M, the latter two conditions both say that 2M/Ri must be 
close to 1. Once ro, ri, i?0) -^^ij and rh(ro) have been chosen, \W-\ has to be 
chosen sufficiently large where we refrain from making this precise here. Any 
initial distribution on {r > ro} which satisfies (|3.6p and (j4.6p is admissible 
in Theorem 14.11 and it is easy to see that there exist such data. 

Since the various parameters which enter the definition of our class of 
admissible data are defined in terms of inequalities, the set of data has "non- 
empty interior" , in the sense that sufficiently small perturbations of initial 
data in the "interior" of this set belong to it as well. 

The crucial step in the proof of the remainder of Theorem 14.11 is to show 
that all particles move towards r = 2M with the stated estimate. As a first 
step note that by ()3.ip and ()3.3p . 



r - 2M 

(// - X)(t, r) > In , r > 2M. 

r 

Together with the control for the radial momentum of the particles this 
implies that along any characteristic in the matter support, 

^ ^ g(/.-A)(s,r)^ < _(J^{^i-\){s,r) < - 2M 

E ~ ~ r 

as long as r > 2M; C > is determined by the initial data parameters. 
Integrating this differential inequality proves the support estimate in the 
theorem. 
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The spherically symmetric, asymptotically flat Einstein- Vlasov system 
has a wide variety of static solutions with finite ADM mass and compact 
matter support, cf. [36 \ I4 H Particularly interesting examples of initial 
data to which our results apply are obtained if the matter for r < tq is 
represented by such a static solution. 

Corollary 4.2 Let fg he a static solution of the spherically symmetric, 
asymptotically flat Einstein- Vlasov system with finite ADM mass Mg > 
and finite radius > 0. Define ro = rg, let ri > rg be arbitrary, M = ri/2, 
and Mout = M — M^; the latter quantity is positive. Then one can construct 
data f on {r > r^} such that the solution of the Einstein- Vlasov system 

o 

launched by fg -\- f has the properties stated in Theorem \4.1\ it exists for all 
t > and r > 0, and it coincides with the static solution fg for all r < 7^(t) 
and t >0. 

It is at this point that the choice 8/9 in (j4.9p is relevant; for the proof 
of Theorem 14.11 any positive constant less than 1 would do. But in |4| it is 
shown that for any steady state of the spherically symmetric Einstein- Vlasov 
system the condition 2m{r)/r < 8/9 holds for all r > 0, and this bound is 
actually sharp, cf. [3j. 

Given the fact that the solutions described in Theorem [4T] undergo gravi- 
tational collapse and form a black hole it may seem strange that in the center 
of such a solution a steady state comfortably sits for all t > 0. But this is 
of course due to the fact that t is Schwarzschild time. In the region r < 2M 
the solution can be extended using a different time coordinate, and if the 
latter is properly chosen then the outer matter will all pass within r < 2M, 
it will crash into the steady state, and all the matter will finally collapse 
into a spacetime singularity. To support all phases of this evolution with 
rigorous theorems is one of our projects for current and future research, cf. 
Section 15.31 

On the other hand, the fact that the steady state sits undisturbed in the 
center as long as none of the outer matter reaches it is easily understood from 
a physics point of view, and just as easily proven. In spherical symmetry 
particle orbits within a certain radius r are not influenced by matter which 
is at strictly larger radii. Hence adding the outer matter shell does not 
change the fact that the steady state satisfies the Einstein- Vlasov system on 
{r < ro} as long as the outer matter remains outside that region. 
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4.3 Weak cosmic censorship and the formation of a black 
hole 

The result of the previous section shows in particular that no particle and 
no light ray can escape from the region {r < 7*(t)} C {r < 2M}. Since 
outside of this domain the solution is global, this means in particular that 
no causal curve originating in a possible singularity can reach the region 
{r > 2M}. However, all this refers to the formulation of the Einstein- 
Vlasov system in Schwarzschild coordinates, and we are not (yet) allowed to 
conclude that our spacetime satisfies the weak cosmic censorship conjecture 
and that a black hole forms. In the present subsection we address these 
questions in a coordinate- independent manner. We start with showing that 
spacetimes as obtained in the previous subsection satisfy the weak cosmic 
censorship conjecture which in heuristic terms says that no singularities can 
be observed from infinity even if the observations are allowed to continue 
indefinitely. Expressed more precisely we have to show: 

Proposition 4.3 The spacetimes obtained in Theorem \4-l\ possess a com- 
plete future null infinity. 

The concept of future null infinity usually refers to a conformal compacti- 
fication of the spacetime under investigation which attaches a boundary at 
infinity, cf. [HH 11.1]. We prefer to follow Christodoulou [20l p. A26] for the 
definition of the term "possess a complete future null infinity". 

The set Bq := {(0,r) | r < has compact closure, and the boundary 
of its causal future is given by the radially outgoing null geodesic 71 
starting at we recall that Ri is the outer radius of the initial matter 
support in Theorem 14. 1[ By that theorem, 71 is future complete. Consider 
now a domain B := {(0,r) | r < R2} with R2 > Ri and the boundary 
C~ of its domain of dependence which is given by the radially ingoing null 
geodesic 72 starting at R2- According to the definition in [20] we must show 
that the affine length of 72 , measured from the intersection of 71 and 72 , goes 
to infinity as R2 ^ 00. In doing so the affine parameterization of 72 must 
be normalized in such a way that its tangent at (0, -R2) equals the vector 
T — N where T is the future directed unit normal to the initial hypersurface 
{t = 0}, i.e., T = (e~^(°'-^2)^ 0,0,0), and N is the outward unit normal to B 
in the initial hypersurface, i.e., N = (0,e~^^^'^^\0,0). 

It turns out that already the affine length of 72 between the intersections 
with 71 and with the line r = Ri goes to infinity as R2 ^ 00. But on the 
region {r > Ri} we have according to Theorem 14.11 vacuum with metric 
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given by 

r 

When parameterized by coordinate time radial null geodesies in the region 
{r > Ri} satisfy the estimates 



1 <1 = e''-^ 

Ri - r 



dr 
lit 



< 1. 



Let (T*,i?*) denote the point where 71 and 72 intersect. By the above 
estimate, 71 (t) < Ri + t and 72(0 > R2 — t which implies that T* > 
{R2 — Ri)/2 and hence 



Consider now an affine parameterization r {t, r, 6, ip){T) of 72, r > 0, with 
(t, r)(0) = (0, i?2)- Since 72 is radial, 9 = 7r/2, if = 0. Since 72 is null, 

and since 72 is ingoing, i.e., dr/dr < 0, we find that 

dt \ ,, dr 
e 

dr dr 

By the geodesic equation, 



0. (|)^ a,. 



dr 



2 



= A7rre'^\^—j [-p - p - 2j] = 0, (4.10) 

i.e., dr/dr = const =: a as long as 72 is in the vacuum region {r > Ri}. 
The normalization condition mentioned above requires that 

(0) = (^_eM0,ii2)-M0,i?2)^^^^ = (^e-M(o,R2)^_g-A(o,ii2)^ 

which means that we should choose a := —e~'^^^'^^\ Let ti and r2 be the 
values of the affine parameter such that (t,r)(ri) is the intersection point of 
72 with 71 and {t, r)(T2) is the intersection point of 72 with the line r = Ri, 
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i.e., 'r(ri) = R* and r(r2) = Ri- For the affine length L of the corresponding 
piece of 72 we find that 

fT'z fRl 
Jti Jr* O- 

> (l- ^\ R2-R1 A _ 2M\ 

- V 7 2 V ^1 y 

— > 00 as i?2 ^ CO- 

This proves the claim of Proposition 14.31 in the sense of [20] . 

We now turn to the question whether Theorem 14 . II also implies the forma- 
tion of a black hole in some appropriate, coordinate-free sense. A maximal 
development of Cauchy data is said to contain a black hole if future null in- 
finity 2^ is complete and the causal past J^{Z^) of future null infinity has 
non-empty complement, cf. [23l Sect. 12]. Intuitively this says that no causal 
curve, i.e., no particle trajectory or light ray, originating in the complement 
of J~{I~^) can reach future null infinity, and that such trapped causal curves 
really exist. In the spacetimes obtained in Theorem 14. II the null geodesic 7* 
does not reach future null infinity, and since we have shown that the latter 
is actually complete, 7* cannot reach future null infinity in any extension 
(such as the maximal development) either. Since is complete for the 
spacetimes we obtained, 7*, in order to reach in the maximal develop- 
ment, would have to enter the outer region D; -i*{T*) D. But following 
7* backwards we see that it must have stayed in D for r < r* and hence 
cannot have reached the region {r > 2M} to begin with. 

In other words, we have in Theorem ST] constructed a globally hyperbolic 
spacetime M which has a complete T'^ and in which the complement of 
J~{X~^) is non-empty. If we consider the maximal development M. of the 
same Cauchy data, it has the same, complete T"*", and since Ai <Z M, the 
complement of J~(I^) in M is non-empty as well. Hence the following is 
true. 



Proposition 4.4 Initial data as specified in Theorem 4-1 lead to the for- 
mation of a black hole in the following sense. The spacetimes M obtained 
in Theorem \4-l\ possess a complete future null infinity I"*" , and its causal 
past J~{I'^) has non-empty complement in the maximal development of the 
data. 
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5 Concluding remarks 



In this final section we want to discuss the question to which extent the 
results from the previous section really depend on the particular matter 
model which we employed, we want to compare these results to previous 
results for other matter models, and we want to discuss future perspectives. 

5.1 General matter models 

The issue of gravitational collapse and in particular the validity of the cosmic 
censorship conjecture should of course be addressed not only for one par- 
ticular matter model like the collisionless gas. Indeed, some key arguments 
in the proof of Theorem 14.11 do depend only on certain general properties 
of the matter model and not on its specific nature. We briefly discuss this 
issue. To this end, let 

p := e-^/^Too, p := e'^^Tu, j := -c-^-^Tqi. (5.1) 

Firstly, we assume that the following two conditions are satisfied. 

• The dominant energy condition holds. (DEC) 

• The radial pressure p is non-negative. (NNP) 

In general relativity the dominant energy condition (DEC) is the main crite- 
rion which a matter model must satisfy in order to be considered realistic, cf. 
[29j . The non- negative pressure condition (NNP) is a standard assumption 
for most astrophysical models. These two criteria imply that 

< p < p and \j\ < p, (5.2) 

cf. [29] and [35] • Furthermore, by (DEC) any world line {s,R{s)) of a ma- 
terial particle or a photon satisfies the estimate 

< g(M-A)(s,fl{s)) 

SO that locally the speed of energy flow is less than or equal to the speed 
of light. In addition we need to assume certain a-priori information on the 
behavior of the solutions of the Einstein- matter equations, namely that for 
solutions launched by data from a suitable class 

• 7^ defined by (j4.2p exists on [0, oo[, and the solution exists on D 

(GLO) 



dR 
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• There exists a constant ci > such that p < —cij in D. (GCC) 

The role of the "global existence condition" (GLO) is obvious. The "grav- 
itational collapse condition" (GCC) is crucial for our method of proof, in 
particular, together with (|5.2|) it implies that j < in D, i.e., the matter is 
ingoing for all times. Notice also that in (j4.8p condition (GCC) was used. 
We emphasize that for Vlasov matter the conditions (DEC) and (NNP) hold 
always, while (GLO) and (GCC) follow via (j4.7p from a suitable restriction 
of the initial data. 

For a general matter model satisfying these assumptions not all the re- 
sults in Theorem 14.11 can be obtained but the following still holds: 

lim 7^*^(5) < 00 and lim fJ.{t,r) = —00 for lim 7^(s) < r < ri 

s — ^00 t — ^00 s — ^00 

for some ri > lims_>oo 7^(s)- If ^* and 7* are defined as in Theorem 14.11 
then 

lim 7*(s) < 00, 

s^oo 

and every radially outgoing null geodesic 7 with 7(0) > r* is future complete 
with lims^oo7('5) = 00. The analysis in Section [4.31 applies to the resulting 
spacetimes as well, in particular since they again are vacuum for r > Ri. 

5.2 Related results 

Due to the inherent difficulties of the Einstein field equations progress to- 
wards understanding the issues of cosmic censorship and the formation of 
black holes has up to now been restricted to the case of spherical symmetry. 
At least this assumption is made in all the papers we mention below. 

The most complete understanding of the issues at hand has so far been 
obtained for the case where the Einstein equations are coupled to a massless 
scalar field, cf . [H HSl [H ttZl [H [H [20] . The final resuh of these investiga- 
tions is that weak and strong cosmic censorship hold for the Einstein-scalar 
field system. A crucial step was the investigation in [16] where explicit con- 
ditions on the initial data were formulated which guarantee the formation 
of a trapped surface. Although the conditions we impose for Theorem 14.11 
are reminiscent of the ones in [16], there is also an important difference. In 
|16j the admissible data cover the full range of 2m/ r g]0, 1[, whereas for our 
result we need that initially 2m /r is close to one in the outer matter ring. 

However, there may be a better reason for this difference than just the 
limited abilities of the present authors. The spherically symmetric Einstein- 
Vlasov system can exhibit at least the following qualitatively different types 
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of behavior. Firstly, for small data the solution disperses in the sense of 
Theorem l3.2[ Secondly, the system has a tremendously rich family of steady 
states, cf. [Ml suns] and also [8] concerning the possible shapes these steady 
states can take. Numerical investigations show clearly that there are both 
stable and unstable steady states, cf. [7]. There is also strong numerical 
evidence [7] that the system has time-periodic or at least almost periodic 
solutions, which is the third type of solution behavior. The fourth and 
last type of solution behavior is gravitational collapse and formation of a 
black hole as shown by Theorem 14.11 To our knowledge only dispersion 
and gravitational collapse are known as possible solution behaviors for the 
Einstein-scalar field system, and this wider range of qualitative solution be- 
haviors for the Vlasov matter model may explain why the condition needed 
to force gravitational collapse is more restrictive than for the scalar field. 

In passing we also note that the primary motivation for coupling the 
Einstein equations to a scalar field is, according to [20j, to capture the wave 
nature of non-symmetric vacuum solutions to the Einstein equations in a 
spherically symmetric set up while still enjoying the simplifications which 
the latter symmetry assumption entails. As mentioned above, the Vlasov 
matter model is actually used in astrophysics to describe galaxies or globular 
clusters. 

Another matter model which features prominently in the history of the 
concepts of gravitational collapse and black hole is a fluid with pressure zero 
which is usually termed dust. The analysis in [32] has definitely shaped the 
overall picture of gravitational collapse to a large extent. The physical reason 
for neglecting pressure is the intuition that once the collapse is sufficiently 
advanced, gravity will dominate all other forces, including pressure. One 
important mathematical advantage of this matter model is that one can use 
coordinates which are co-moving with the matter. On the other hand, dust 
produces shell-crossing singularities which hinder the analysis of the real 
issues, and dust can lead to naked singularities as shown in [13]. Whether 
cosmic censorship will eventually be established for Vlasov matter or not, 
solutions of the Einstein- Vlasov system do have non-trivial pressure, and 
the system has so far not been shown to produce naked singularities. 

In making the last statement we are fully aware of [50]. There it is 
claimed that numerical simulations of the axially symmetric Einstein- Vlasov 
system can lead to naked singularities. However, there is reason to believe 
that the matter model which was actually simulated in [50] was dust and 
not Vlasov, cf. [i5] . 

To conclude this comparison with previously known results we mention 
that in [45] a continuity argument was used to show that there exist initial 
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data for the spherically symmetric Einstein- Vlasov system which lead to the 
formation of trapped surfaces. Combining this with the results of [21j and 
[23) implies that there exist data which lead to the formation of a black 
hole. Due to the method of proof in 05] this analysis does not give explicit 
conditions on the data which guarantee this type of behavior. Our approach 
does produce such explicit conditions on the data, and these data are stable 
against perturbations which are small in a suitable sense. 

5.3 Open problems and future perspectives 

To conclude we mention some open problems and possible perspectives for 
future research. An immediate question is whether the geodesic 7* in The- 
orem 14.11 is future complete; the question whether an event horizon has 
complete generators is of general interest. This issue is closely related to 
the question whether in the limit t — > 00 all the matter ends up in the 
region {r < 2M}, i.e., whether lim.t^oon^{t,2M) = 2M. It is clear from 
the geodesic equation for 7* that this geodesic can not be complete if it is 
running in vacuum all the time, cf. (14.101) . which means that at least some 
matter must cross r = 2M if 7* should be complete. These questions are 
currently under investigations by two of the authors, and we believe that for 
suitably restricted data indeed limt^oo rn{t, 2M) = 2M. 

The next logical step will be to analyze the situation corresponding to 
Theorem l4.1l in a coordinate system which has the potential to cover regions 
of spacetime containing trapped surfaces and which may reach all the way 
to the spacetime singularity which forms at the center. Preliminary steps in 
this direction look promising. 

Much more demanding is the question of how to relax the conditions in 
Theorem 14.11 so that the numerical observations reported in [7] where the 
perturbation of unstable steady states leads to the formation of black holes 
are covered. An analytic understanding of the stability properties which 
were observed numerically is a further open problem. An answer to this 
could also help to explain the phenomenon of critical collapse, cf. [TJ EU HI] 
for relevant numerical results for the Einstein- Vlasov system and [28j for a 
general discussion of this issue. 

As pointed out above the question whether weak cosmic censorship holds 
for the spherically symmetric Einstein- Vlasov system is related to the ques- 
tion whether this system has global solutions in Schwarzschild coordinates 
or not. Useful estimates in the direction of global existence in Schwarzschild 
coordinates, which go beyond what was reported in Section [3l were estab- 
lished in ^ , but the problem remains open. And when thinking about weak 
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cosmic censorship one should definitely keep the possibility in mind that 
global existence in Schwarzschild coordinates could be violated for an initial 
data set "of measure zero" , while weak cosmic censorship could still be true 
for the spherically symmetric Einstein- Vlasov system. 

Eventually one will wish to go beyond spherical symmetry, and any ex- 
tension of the results mentioned in these notes to for example the case of 
axial symmetry is in our opinion a challenging and worthwhile problem. 
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